High energy cosmic-rays from gamma-ray burst sources: A stronger case 

Eli Waxman 1 

ABSTRACT 

The suggested association between the sources of 7-ray bursts (GRB's) and the 
sources of ultra-high energy cosmic rays (UHECR's) is based on two arguments: (i) 
The constraints that UHECR sources must satisfy to allow proton acceleration to > 
10 20 eV are similar to those inferred for GRB sources from 7-ray observations, and 
(ii) The average energy generation rate of UHECR's is similar to the 7-ray generation 
rate of GRB's. We show that recent GRB and UHECR observations strengthen both 
arguments, and hence strengthen the suggested association. 

Subject headings: acceleration of particles — cosmic-rays — gamma-rays: bursts — gamma 
rays:theory 

1. Introduction 

The widely accepted interpretation of the phenomenology of 7-ray bursts (GRB's), bursts of 
0.1 MeV-1 MeV photons lasting for a few seconds (see Fishman & Meegan 1995, for review), is that 
the observable effects are due to the dissipation of the kinetic energy of a cosmologically distant, 
relativistically expanding wind, a "fireball," whose primal cause is not yet known (for reviews see 
Piran 2000; Meszaros 2002; Waxman 2003a). Waxman (1995a), Milgrom & Usov (1995) and 
Vietri (1995) have suggested that ultra-high energy, > 10 19 eV, cosmic rays (UHECR's) may be 
produced in GRB sources. The model suggested in (Waxman 1995a) was based on two arguments. 
First, it was shown that the constraints imposed on the relativistic wind by the requirement that it 
produces observed GRB characteristics are similar to the constraints imposed on such a wind by the 
requirement that it would allow proton acceleration to > 10 20 eV. Second, the energy generation 
rate of 7-rays by GRB's was shown to be similar to the energy generation rate required to account 
for the observed UHECR flux (Waxman 1995a, b). 

The origin of UHECR's is one of the most exciting open questions of high energy astrophysics 
(Bhattacharjee & Sigl 2000; Nagano & Watson 2000). The extreme energy of the highest energy 
events poses a challenge to models of particle acceleration. Since very few known astrophysical 
objects have characteristics indicating that they may allow acceleration of particles to the observed 
high energies (Waxman 2003b), the question of whether GRB's are possible UHECR sources is of 
great interest. Moreover, since the GRB model for UHECR production makes unique predictions, 
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which differ from those of other models (see Waxman 2001, and discussion in § 4), the design and 
analysis of future large area UHECR experiments may be affected by the answer to this question. 

The detection over the past few years of "afterglows," delayed low energy (X-ray to radio) 
emission of GRBs (for review see Kulkarni et al. 2000) confirmed the cosmological origin of 
the GRB's, through redshift determination of GRB host-galaxies, and confirmed standard model 
predictions of afterglows that result from the collision of an expanding fireball with its surrounding 
medium (e.g. Piran 2000; Meszaros 2002; Waxman 2003a). In addition to providing strong support 
to the fireball model, these observations also provide new constraints on fireball model parameters 
and more accurate information on the redshift distribution of GRB sources. Recently, new data 
on the spectrum and flux of UHECR's was presented by the HiRes experiment (Abu-Zayyad et al. 
2002), providing improved constraints on the generation rate and spectrum of UHECR's (Bahcall 
& Waxman 2003). Here we discuss the implications of these new GRB and UHECR observations 
to the GRB model for UHECR production. 

In § 2.1 we briefly describe the model proposed in (Waxman 1995a) for proton acceleration in 
GRB fireballs. The main goal of this section is to identify the key constraints that the relativistic 
wind parameters need to satisfy, Eqs. (1) and (2), in order to allow acceleration of protons to 
> 10 20 eV (A more detailed and pedagogical description of the model is given in Waxman 2001). 
The association of GRB and UHECR sources was motivated mainly by the fact that these con- 
straints were similar to those inferred, based on independent physical arguments, from > 1 MeV 
7-ray observations. Recent lower energy afterglow observations provide new constraints on model 
parameters, which are independent of the 7-ray constraints. In § 2.2 we show that the afterglow 
constraints are similar to the 7-ray constraints, and in fact imply parameter values which are more 
favorable for the acceleration of protons to > 10 20 eV. In § 2.3 we compare our results to those 
of other authors. In particular, we show that results derived by several authors, arguing that the 
maximum proton energy is too low (Gallant & Achterberg 1999; Achterberg et al. 2001) or that 
the proton spectrum is too steep (Ostrowski & Bednarz 2002) to account for the observed UHECR 
spectrum, are not applicable to the model proposed in (Waxman 1995a). 

The determination of GRB redshifts, made possible by afterglow detection, lead to significant 
changes in, and to significantly reduced uncertainties of, the estimates of both GRB rate and 
average 7-ray energy release (per single GRB). We show in § 3.1 that the local energy production 
rate in 7-rays by GRB's, inferred using recent redshift measurements, is similar to the pre-afterglow 
estimate. In § 3.2 we show that the local rate of energy production in UHECR's inferred using 
recent UHECR observations is consistent with, although more accurate than, earlier estimates 
(Waxman 1995b), and hence comparable to the energy production rate in 7-rays by GRB's. In 
§ 3.3 we demonstrate that recent claims to the contrary (Stecker 2000; Berezinsky et al. 2002), in 
particular that the UHECR energy production rate is 3 order of magnitude higher than the GRB 
7-ray production rate (Berezinsky et al. 2002), are erroneous. 

Our main results are summarized in § 4. 
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2. Proton acceleration in GRB fireballs 

2.1. Brief description of the model 

General phenomenological considerations, based on 7-ray observations, indicate that, regard- 
less of the nature of the underlying sources, GRB's are produced by the dissipation of the kinetic 
energy of a relativistic expanding fireball. A compact source, r$ ~ 10 7 cm, produces a wind, char- 
acterized by an average luminosity L ~ 10 52 erg s" 1 and mass loss rate M. At small radius, the 
wind bulk Lorentz factor, T, grows linearly with radius, until most of the wind energy is converted 
to kinetic energy and T saturates at T ~ L/Mc 2 ~ 300. Variability of the source on a time scale 
At ~ 10 ms, resulting in fluctuations in the wind bulk Lorentz factor T on a similar time scale, 
results in internal shocks in the ejecta at a radius r ~ ~ T 2 cAt 3> r^. It is assumed that internal 
shocks reconvert a substantial part of the kinetic energy to internal energy, which is then radiated 
as 7-rays by synchrotron and inverse-Compton radiation of shock-accelerated electrons. At a later 
stage, the shock wave driven into the surrounding medium by the expanding fireball ejecta leads 
to the emission of the lower-energy afterglow. 

The observed radiation is produced, both during the GRB and the afterglow, by synchrotron 
emission of shock accelerated electrons. In the region where electrons are accelerated, protons are 
also expected to be shock accelerated. This is similar to what is thought to occur in supernovae 
remnant shocks, where synchrotron radiation of accelerated electrons is the likely source of non- 
thermal X-rays, and where shock acceleration of protons is believed to produce cosmic rays with 
energy extending to ~ 10 15 eV (e.g. Koyama et al. 1995; Berezhko et al. 2002; Volk et al. 2002). 
Thus, it is likely that protons, as well as electrons, are accelerated to high energy within GRB 
fireballs. 

The internal shocks within the expanding wind are expected to be mildly relativistic in the 
wind rest frame, due to the fact that the allowed range of Lorentz factor fluctuations within the wind 
is from few x 10 2 (the lower limit required to avoid large optical depth) to few x 10 3 (the maximum 
Lorentz factor to which shell acceleration by radiation pressure is possible, e.g. Waxman 2003a). 
This implies that the Lorentz factors associated with the relative velocities are not very large. 
Since internal shocks are mildly relativistic, we expect our understanding of non-relativistic shock 
acceleration to apply to the acceleration of protons in these shocks. In particular, the predicted 
energy distribution of accelerated protons is expected to be dn p /dE p oc E p 2 (Axford et al. 1977; 
Bell 1978; Blandford & Ostriker 1978), similar to the predicted electron energy spectrum, which 
is consistent with the observed photon spectrum. 

Several constraints must be satisfied by wind parameters in order to allow proton acceleration 
to high energy E p . We summarize below these constraints. The reader is referred to (Waxman 
1995a, 2001) for a detailed derivation. The requirement that the acceleration time be smaller than 
the wind expansion time (which also implies that the proton is confined to the acceleration region 
over the required time) sets a lower limit to the strength of the wind magnetic field. This may 
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be expressed as a lower limit to the ratio of magnetic field to electron energy density (Waxman 
1995a), 

u B /u e > 0.02r| 5 £2 2oL ~i 2; (1) 

where E p = 10 20 £ , Pi 20 eV, T = 10 2 - 5 r 2 .5 and L 7 = 10 52 L 7i 52erg/s is the wind 7-ray luminosity A 
second constraint is imposed by the requirement that the proton acceleration time be smaller than 
the proton energy loss time, which is dominated by synchrotron emission. This sets an upper limit 
to the magnetic field strength, which in turn sets a lower limit to T (Waxman 1995a; Rachen & 
Meszaros 1998) 

r > l30E 3 p %Arl /4 . (2) 

Here, At = 10 _2 At_2 s. As explained in (Waxman 1995a), the constraints Eq. (1) and Eq. (2) 
hold regardless of whether the fireball is a sphere or a narrow jet (as long as the jet opening angle is 
> 1/r). The luminosity in Eq. 1 is the "isotropic equivalent luminosity", i.e. the luminosity under 
the assumption of isotropic emission. 

Internal shocks within the wind take place at a radius ~ T 2 cAt. The constraint of Eq. (1) is 
independent of At, i.e. independent of the internal collision radius, while the constraint of Eq. (2) 
sets a lower limit to the collision radius for a given At. This implies that protons may be accelerated 
to > 10 20 eV regardless of the value of At, which may range from the dynamical time of the source 
(At ~ 1 ms) to the wind duration (At ~ 1 s), provided the magnetization and Lorentz factor are 
sufficiently large, following Eqs. 1 and 2. 

At large radii the external medium affects fireball evolution, and a "reverse shock" is driven 
backward into the fireball ejecta and decelerates it. For typical GRB fireball parameters this shock 
is also mildly relativistic (e.g. Waxman 2003a), and its parameters are similar to those of an 
internal shock with At ~ 10 s. Protons may therefore be accelerated to > 10 20 eV not only in the 
internal wind shocks, but also in the reverse shock (Waxman & Bahcall 2000; Waxman 2001). 
This implies that proton acceleration to > 10 20 eV is possible, provided the constraints of Eqs. 1 
and 2 (with At ~ 10 s) are satisfied, also in (the currently less favorable) scenario where GRB 
7-rays are produced in the shock driven by the fireball into the surrounding gas, rather than by 
internal collisions (as suggested, e.g., in Dermer 1999). 

The constraints given by Eqs. (1) and ( 2) are remarkably similar to those inferred from 
7-ray observations, based on independent physical arguments: T > 300 is implied by the 7-ray 
spectrum by the requirement to avoid high pair-production optical depth, and magnetic field close 
to equipartition, us/ue ~ 0.1, is required in order to account for the observed 7-ray emission 
(Piran 2000; Meszaros 2002; Waxman 2003a). This was the basis for the association of GRB's 
and UHECR's suggested in (Waxman 1995a). In the following sub-section we discuss the new 
constraints on model parameters implied by afterglow observations, and their implications. 
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2.2. Implications of afterglow observations 

Afterglow observations lead to the confirmation of the cosmological origin of GRBs and con- 
firmed standard model predictions of afterglow that results from synchrotron emission of electrons 
accelerated to high energy in the highly relativistic shock driven by the fireball into its surrounding 
gas. Afterglow observations provide therefore strong support for the underlying fireball scenario. In 
addition, afterglow observations provide important information on the values of model parameters 
that enter the constraints given by Eqs. (1) and (2). 

Prior to the detection of afterglows, it was commonly assumed that the farthest observed 
GRB's lie at redshift z ~ 1 (Mao & Paczyhski 1992; Piran 1992). Based on afterglow redshift 
determinations, we now know that detected GRB's typically lie at farther distances (e.g. Bloom, 
Frail & Kulkarni 2003). This implies that the characteristic GRB luminosity is higher by an order 
of magnitude compared to pre-afterglow estimates, L 7 ss 10 52 erg/s instead of L 7 ss 10 51 erg/s. This 
relaxes the constraint on magnetic field energy fraction given by Eq. (1). The implications of the 
revised GRB redshift distribution to the inferred GRB energy production rate are discussed in 
detail in § 3. 

In several cases, fast follow up afterglow observations allowed the detection of radio and optical 
emission from the reverse shock (Zhang, Kobayashi & Meszaros 2003; Soderberg & Ramirez- 
Ruiz 2003, and references therein). These observations provide direct information on the plasma 
conditions in the reverse shock, where acceleration of protons to high energy may take place (see 
§ 2.1). Two major conclusions were drawn from the analysis of the early optical and radio reverse 
shock emission. First, lower limits to the initial fireball Lorentz factors were inferred, in the range 
of T > 100 to T > 1000 (Zhang, Kobayashi & Meszaros 2003; Soderberg & Ramirez-Ruiz 2003). 
Second, the magnetic field in the reverse shock was inferred to be close to equipartition, that is 
UB/u e was inferred to be of order unity (Waxman & Draine 2000; Zhang, Kobayashi & Meszaros 
2003). Early afterglow observations provide therefore constraints on T and on UB/u e which arc 
(i) Independent of the constraints derived from 7-ray observations; (ii) Consistent with the 7-ray 
constraints; and (hi) Are remarkably similar to the constraints of Eqs. (1) and (2), that need to be 
satisfied in order to allow proton acceleration to > 10 20 eV. 

2.3. Comparison with other authors 

Gallant & Achterberg (1999) and, more recently, Achterberg et al. (2001) have considered 
particle acceleration by the ultra-relativistic, T ~ 300 shock driven by the fireball into its sur- 
rounding medium. They argue that protons can not be accelerated in this external ultra-relativistic 
shock to ultra-high energy. Regardless of whether or not this claim is valid, it is irrelevant for 
the model proposed in (Waxman 1995a) and discussed in § 2.1, where proton acceleration takes 
place in internal (reverse) mildly-relativistic shocks. Similarly, the claims in Ostrowski & Bednarz 
(2002), that the spectrum of protons accelerated in ultra-relativistic shocks is much steeper than 
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dnp/dE p oc E p 2 , the spectrum expected for sub-relativistic shocks and required to account for the 
observed UHECR spectrum (see § 3.2), are not applicable to the model proposed in (Waxman 
1995a). 

It should be emphasized here, that non-relativistic collisionless shocks are observed in many 
types of astrophysical systems, and that the theoretical understanding of particle acceleration in 
such shocks (Drury 1983; Blandford & Eichler 1987) is more developed than in the case of 
relativistic collisionless shocks. The GRB model for particle acceleration, described in § 2.1, relies 
on our understanding of acceleration in non-relativistic collisionless shocks, and therefore is not 
subject to the uncertainties described in the preceding paragraph, which are related to acceleration 
in ultra-relativistic shocks. 

Association of GRB's and UHECR's has been also suggested by Vietri (1995) and by Milgrom 
&; Usov (1995), who noted that GRB's may accelerate protons to > 10 20 eV energy. While Milgrom 
& Usov did not suggest an acceleration mechanism, the criticism of Gallant et al. (1999) may be 
relevant to the mechanism proposed by Vietri. Their main point is that the fractional energy gain 
per shock crossing is of order unity for highly relativistic shocks, rather than of order T 2 as suggested 
in Vietri (1995). This implies that the acceleration time in ultra-relativistic shocks is much longer 
than estimated by Vietri (1995). However, this conclusion does not necessarily imply that the 
acceleration process is ineffective. In their estimate of the acceleration time Gallant & Achterberg 
(1999) and Achterberg et al. (2001) used an up-stream magnetic field amplitude of 1/iG, typical 
to the inter-stellar medium. GRB observations imply that this pre-shock magnetic field must be 
amplified by many orders of magnitude in the GRB shock (Gruzinov & Waxman 1999), and it 
is therefore far from clear that the magnetic field value relevant for particle deflection up-stream 
is the un-perturbed pre-shock field. The up-stream magnetic field may be amplified ahead of the 
shock by, e.g., the streaming of high energy particles (Bell & Lucek 2001; Dermer 2002), in which 
case acceleration to ultra-high energy is possible. 



3. UHECR energy generation rate and spectrum 

3.1. The GRB energy generation rate 

The GRB model for UHE cosmic-ray production was suggested prior to the detection of after- 
glows. Estimates of the rate of GRB's were based at that time on the 7-ray flux distribution, and 
ranged from ~ 3/Gpc 3 yr (Piran 1992) to ~ 30/Gpc 3 yr (Mao & Paczyhski 1992). The estimated 
average 7-ray energy release in a single GRB, based on a characteristic peak flux of ~ 10 51 erg/s 
(e.g. Fishman h Meegan 1995), was ~ 10 52 erg. These estimates were subject to large uncertain- 
ties, since the 7-ray luminosity function as well as the evolution of GRB rate with redshift were 
poorly constrained. Based on the rate and energy estimates, the rate of 7-ray energy generation 
by GRB's was estimated to be ~ 10 44 erg/Mpc 3 yr. The determination of GRB redshifts, which was 
made possible by the detection of afterglows, allows a more reliable estimate. 
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Most of the GRB's are observed from z > 1, since they can be detected out to large redshift. 
This implies that the GRB rate density at z > 1 is better constrained by the observations than 
the local, z = 0, rate. The inferred local rate depends on the assumed redshift evolution. It is now 
commonly believed that the GRB rate evolves with redshift following the star-formation rate, based 
on the association of GRB's with type Ib/c supernovae. This association, which was originally 
motivated by the temporal and angular coincidence of GRB980425 and SN1998bw (Galama et 
al. 1998), has gained significant support from the identification of a SN1998bw-like spectrum 
in the optical afterglow of GRB030329 (Stanek et al. 2003; Hjorth et al. 2003). It is also 
supported by evidence for optical supernovae emission in several GRB afterglows (Bloom 2003). 
Adopting the assumption, that the GRB rate follows the redshift evolution of the star formation 
rate, the local (z = 0) GRB rate density was inferred by Schmidt (2001) to be Rgrb(z = 0) rs 
0.5 x 10~ 9 Mpc~ 3 yr _1 . A similar result was later obtained, under similar assumptions, by Perna, 
Sari & Frail (2003). 

The local rate density determined by Schmidt (2001) and by Perna, Sari & Frail (2003) is 
uncertain due to uncertainties in the determination of the evolution of the star formation rate: 
A faster evolution (rate increase) with redshift implies a lower local, z = 0, GRB rate. In their 
analysis, both Schmidt (2001) and Perna, Sari & Frail (2003) have not used, however, the detailed 
information provided by BATSE on the observed distribution of GRB peak fluxes, and the detailed 
shape of the observed GRB redshifts distribution (Schmidt 2001, has used only the value of 
< V/V max >). In a more detailed analysis, taking into account these constraints, Guetta, Piran & 
Waxman (2003) have shown that observations allow to discriminate between different assumptions 
regarding the evolution with redshift of the GRB rate density. Redshift evolution following the 
Rowan-Robinson (Rowan-Robinson 1999) star formation rate evolution (i?GRB io°- 75z U p to 
z = 1) was found consistent with observations, while a more rapid evolution (as assumed, e.g. 
in Porciani & Madau 2001), was found inconsistent. For the Rowan-Robinson evolution, Guetta, 
Piran & Waxman (2003) find Rgrb(z = 0) « 0.5 x 10 _9 Mpc~ 3 yr _1 . Given the current (systematic 
uncertainties in the redshift) data, this rate is accurate to within a factor of a few (Guetta, Piran 
& Waxman 2003). 

The local energy generation rate in 7-rays by GRB's, e 7 , is given by the product of R G rb(z = 0) 
and the average 7-ray energy release in a single GRB, e 7 . Bloom, Frail & Kulkarni (2003) provide 
e 7 for 27 bursts with known redshifts, in a standard rest-frame bandpass, 0.02 MeV to 2 MeV. The 
average is e 7 = 2.9 x 10 53 erg, with estimated uncertainty, due to the correction to a fixed rest-frame 
bandpass, of ~ 20% for individual bursts (and much smaller for the average). In calculating e 7 
from this value of e 7 , the following point should be taken into account. e 7 is the average energy 
for bursts with known redshift, most of which were localized by the BeppoSAX satellite. Since 
BeppoSAX has a higher detection flux threshold than BATSE (see Band 2003; Guetta, Piran & 
Waxman 2003), it is sensitive to ~ 70% of the bursts detectable by BATSE, for which the GRB rate 
Rgrb(z = 0) was inferred. Thus, the energy generation rate by bursts detectable by BeppoSAX is 

^[ao2MeV,2MeV] ~ 0.7R GRB (z = 0)e 7 = 10 44 erg Mpc" 3 yr" 1 . (3) 
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As mentioned above, the main uncertainty in determining e 7 is related to the uncertainty in the 
local GRB rate, due to which the value given in Eq. 3 is accurate to within a factor of few. 

The energy generation rate given by Eq. 3 does not reflect the total energy emitted by GRB's 
in 7-rays, rather the energy emitted in the [0.02MeV, 2MeV] band. Emission at higher energy 
has been detected in many GRB's (see, e.g. Gonzalez et al. 2003), implying that the total 7-ray 
energy emission is higher than that limited to the [0.02MeV, 2MeV] band. A bandpass independent 
results may be obtained as follows. GRB 7-ray spectra are well described by broken power-laws, 
dn 1 /dE 1 oc E~ a for £ 7 < E b and E~ P for £ 7 > E b (Band 1993; Preece et al. 1998). The 
observed (redshifted) break energy E b is typically a few hundered keV, j3 ~ 2 and a « 1. The 
observed 7-rays are produced in the fireball model by synchrotron and inverse-Compton emission 
of electrons accelerated to high energy by collisionless shocks in the expanding fireball wind (see 
§ 2 for more detail). The high energy part of the GRB spectrum is produced by energy loss of 
relativistic electrons, accelerated to a power-law distribution, dn e /dE e oc E~ 2 . Taking into account 
the fact that the photon energy is proportional to the square of the electron energy (and hence 
that log(.E 7 ) spans twice the range spanned by log(£' e )), and that the rest frame break energy is 
higher than observed by a factor 1 + z, the rate of energy generation in relativistic electrons is 

El^- - ^ 2M eV,2MeV] « 10 44 erg Mpc" 3 yr" 1 . (4) 

Finally, the following point should be mentioned. The numbers quoted above for the GRB 
rate density and 7-ray energy release are based on the assumption that GRB 7-ray emission is 
isotropic. If, as now commonly believed (see Frail et al. 2001, and references therin), the emission 
is confined to a solid angle Af2 < 4ir, then the GRB rate is increased by a factor (AQ/47r) _1 and 
the GRB energy is decreased by the same factor. However, their product, the energy generation 
rate, is independent of the solid angle of emission. 



3.2. Comparison with UHECR observations 

The cosmic-ray spectrum flattens at ~ 10 19 eV (Bird et al. 1995; Hayashida et al. 1999). 
There are indications that the spectral change is correlated with a change in composition, from 
heavy to light nuclei (Bird et al. 1995; Dawson et al. 1998; Abu-Zayyad et al. 2001). These 
characteristics, which are supported by analysis of Fly's Eye, AGASA and HiRes-MIA data, and 
for which some evidence existed in previous experiments (Watson 1991), suggest that the cosmic 
ray flux is dominated at energies < 10 19 eV by a Galactic component of heavy nuclei, and at UHE 
by an extra-Galactic source of protons. Also, both the AGASA and Fly's Eye experiments report 
an enhancement of the cosmic-ray flux near the Galactic disk at energies < 10 18 - 5 eV, but not at 
higher energies (Bird et al. 1998; Hayashida et al. 1999a). Fly's Eye stereo spectrum is well fitted 
in the energy range 10 17,6 eV to 10 19,6 eV by a sum of two power laws: A steeper component, with 
differential number spectrum J oc E~ 3 - 50 , dominating at lower energy, and a shallower component, 
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J oc E~ 2S1 , dominating at higher energy, E > 10 19 eV. The data are consistent with the steeper 
component being composed of heavy nuclei primaries, and the lighter one being composed of proton 
primaries. 

Bahcall & Waxman (2003) have shown that the observed UHECR flux and spectrum may be 
accounted for by a two component, Galactic + extra-Galactic model. For the Galactic component, 
this model adopts the Fly's Eye fit, 

The spectrum and energy generation rate of extra-Galactic protons are given in this model by 



Ep—j^- = 0.65 x 10 44 erg Mpc" 3 yr" 1 . (6) 

LIU/ p 



The spectral index, 2, is that expected for acceleration in sub-relativistic collisionless shocks in 
general, and in particular for the GRB model discussed in § 2.1. The energy generation rate 
integrated over the energy range of 10 19 eV to 10 21 eV is 



e 



- ru - 3 x 10 44 erg Mpc -3 yr" 1 . (7) 



[10 19 ,10 21 ]eV 



Uncertainties in the absolute energy calibration of the experiments lead to uncertainty of ~ 20% 
in this rate (Bahcall & Waxman 2003). 

The model used in (Bahcall & Waxman 2003) is similar to that proposed in Waxman (1995b). 
The improved constraints on UHECR spectrum and flux provided by the recent observations of 
HiRes do not change the estimates given in (Waxman 1995b) for the energy generation rate and 
spectrum, Eq. (6), but reduce the uncertainties. Comparing Eqs. (6), or (7), and (4) (which, as 
explained in § 3.1, is accurate to within a factor of a few) we find that the rate at which GRB's 
produce energy in accelerated high energy electrons is comparable to the rate at which energy 
should be produced in high energy protons in order to account for the observed UHECR flux. 

It is important to emphasize the following point. As explained in § 3.1, the 7-ray energy 
production rate reflects the rate at which energy is produced by GRB's in relativistic electrons. 
The ratio between the energy carried by relativistic electrons and protons is not known from 
basic principles (observations suggest that in sub-relativistic shocks protons carry ~ 10 times more 
energy than electrons, e.g. Blandford & Eichler 1987; Berezhko et al. 2002). Thus, an exact match 
between the 7-ray and UHECR generation rates should not in general be expected. Rather, the 
two rates are expected to be similar to within an order of magnitude. 



3.3. Comparison with other authors 



Recent claims, that the GRB energy generation rate is too small to account for the observed 
UHE cosmic-ray flux, (Stecker 2000; Berezinsky et al. 2002), are not valid. 
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In (Stecker 2000) it is argued that at most "10% of the cosmic rays observed above 10 20 eV 
can be accounted for by GRB's." The UHECR energy generation rate estimated in (Stecker 2000) 
is similar to the rate given by Eq. (6), as derived in (Waxman 1995b; Bahcall &; Waxman 2003), 
but the GRB 7-ray energy generation rate estimated in (Stecker 2000) is smaller by approximately 
an order of magnitude compared to the rate derived here, Eq. (3). The 7-ray energy generation 
rate by GRB's is underestimated in (Stecker 2000) by an order of magnitude, since it is based on 
an earlier, less accurate estimate of the GRB rate per unit volume, compare (Schmidt 1999) to 
(Schmidt 2001; Frail et al. 2001; Guetta, Piran & Waxman 2003), and neglects GRB 7-rays in 
the 1 MeV to 2 MeV band. 

Berezinsky (2002) claims that the required UHECR generation rate is three orders of magnitude 
larger than the GRB 7-ray energy generation rate. He argues that the energy generation rate of 
UHECR's implied by observations is ~ 3 x 10 46 erg Mpc -3 yr _1 . The discrepancy with our derived 
rate, Eq. (7), arises because Berezinsky assumes that extra-Galactic protons dominate the observed 
flux in the range of 10 17 eV to 10 18 eV. However, most observers attribute the flux in this energy 
range to Galactic cosmic rays. Thus, in our model, Eqs. (5-6), extra-Galactic sources dominate 
the flux only above 10 19 eV, and the fraction of the observed cosmic-ray flux that is contributed 
by the extra-Galactic component at the energy range of 10 eV to 10 18 eV is 



Observations strongly suggest that cosmic-rays above and below 10 eV originate from different 
sources: Galactic sources of heavy nuclei are likely to dominate below 10 19 eV, while extra-Galactic 
proton sources are likely to dominate at higher energy. Thus, by assuming that extra-Galactic 
particles dominate the flux down to 10 17 eV, Berezinsky greatly overestimates the required en- 
ergy for extra-galactic sources. In addition, the GRB energy generation rate is underestimated in 
(Berezinsky et al. 2002) in the same way as in (Stecker 2000). 



The main constraints that a relativistic wind (fireball) need to satisfy to allow proton accel- 
eration to > 10 20 eV are given by Eqs. (1) and (2): The magnetic field energy density ub should 
exceed a few percent of the relativistic electron energy density u e , and the wind Lorentz factor 
r should exceed ~ 10 2 . The similarity of these constraints and the constraints imposed on wind 
parameters, based on independent physical considerations, by 7-ray observations were the basis for 
the association of GRB and UHECR sources suggested in (Waxman 1995a). Afterglow observa- 
tions have shown that the characteristic GRB luminosity is higher than estimated based on 7-ray 
observations alone, 10 52 erg/s instead of 10 51 erg/s, relaxing the constraint of Eq. (1) on ub/u c . In 
addition, early optical and radio afterglow observations provide new constraints on wind parame- 
ters: They imply large Lorentz factors, V > 10 2 to V > 10 3 , and large magnetic field energy density 




(8) 
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in the fireball plasma, UB/u e of order unity (see § 2.2). These constraints are consistent with those 
previously inferred from 7-ray observations, and with the constraints imposed by the requirement 
to allow proton acceleration to > 10 20 eV. 

The local, z = 0, 7-ray energy generation rate by GRB's is ~ 10 44 erg/Mpc 3 yr (see Eqs. (3), 
(4)). Afterglow observations allow to determine this rate to within a factor of a few. The main 
uncertainty is due to uncertainties in the redshift evolution of the GRB rate (see § 3.1). The local 
rate of energy generation in high energy protons, required to account for the observed UHECR flux, 
is determined by observations with smaller uncertainty, and is also « 10 44 erg/Mpc 3 yr (see § 3.2, 
eqs. (6), (7)). 

Prior to direct GRB redshift measurements, which became possible with the detection of 
afterglows, estimates of the GRB rate ranged from ~ 3/Gpc 3 yr (Piran 1992) to ~ 30/Gpc 3 yr 
(Mao & Paczyfiski 1992). These estimates were subject to large uncertainties, since the 7-ray 
luminosity function as well as the evolution of GRB rate with redshift were poorly constrained. The 
measurements of GRB redshifts allow a more accurate estimate, 0.5/Gpc 3 yr at z = 0. Although this 
rate is significantly lower than the pre-afterglow estimates, the estimated rate of energy generation 
in 7-rays by GRBs is similar to the pre-afterglow estimate, which was ~ 10 44 erg/Mpc 3 yr. This is 
due to the fact that the average 7-ray energy release in a single GRB is larger than the pre-afterglow 
estimates, 10 53 - 5 erg instead of ~ 10 52 erg. 

The numbers quoted above for the GRB rate density and 7-ray energy release are based on 
the assumption that GRB 7-ray emission is isotropic. If, as now commonly believed (see Frail et al. 
2001, and references therin), the emission is confined to a solid angle AO, < 4ir, then the GRB rate 
is increased by a factor and the GRB energy is decreased by the same factor. However, 

their product, the energy generation rate, is independent of the solid angle of emission. 

The local GRB rate implies that the rate of GRB's out to a distance from which most protons of 
energy exceeding 10 20 eV originate, ~ 90 Mpc (see figure 2 in Waxman 1995b), is ~ 10~ 3 /yr. The 
number of GRB's contributing to the observed flux at any given time is given by the product of this 
rate and the spread in arrival time of protons, due to the combined effect of stochastic propagation 
energy loss and deflection by magnetic fields (Waxman 1995a). This time spread may be as large as 
10 7 yr for 10 20 eV originating at 90 Mpc distance, implying that the number of GRBs contributing 
to the > 10 20 eV flux at any given time may reach ~ 10 4 (Waxman 2001). The upper limit on the 
strength of the inter-galactic magnetic field, combined with the low local rate of GRB's, leads to 
unique predictions of the GRB model for UHECR production (Miralda-Escude & Waxman 1996; 
Waxman & Miralda-Escude 1996), that may be tested with operating (Abu-Zayyad et al. 2002), 
under-construction (Cronin 1992) and planned (Teshima et al. 1992) large area UHECR detectors. 
In particular, a critical energy is predicted to exist, 10 20 eV < E c < 4 x 10 20 eV, above which a few 
sources produce most of the UHECR flux, and the observed spectra of these sources is predicted to 
be narrow, AE/E ~ 1: The bright sources at high energy should be absent in UHECRs of much 
lower energy, since particles take longer to arrive the lower their energy. The model also predicts 
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the emission of high-energy, > 1 TeV neutrinos (Waxman & Bahcall 1997), a prediction that 
may be tested with operating (AMANDA), under-construction (ANTARES,IceCube, NESTOR) 
and planned (NEMO) large volume neutrino telescopes (see Halzen & Hooper 2002 for review). 
For more detailed discussion of model predictions, see (Waxman 2001). 

At energies > 10 20 eV, the predicted number, N, of events in conventional models is uncertain 
due to the unknown clustering scale, ro, of the sources, cr(iVp re< ji C ted) /-^predicted = 0.9(ro/10 Mpc) 0,9 
(Bahcall & Waxman 2000). For GRB's in particular, the flux above 3 x 10 20 eV is expected to be 
dominated by few sources (Miralda-Escude & Waxman 1996), and hence large deviations from a 
homogeneous source distribution may be expected (see Waxman (2001) for detailed discussion). 
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